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Abstract 

If neutrinos have mass, they can also decay. Decay hfetimes of cosmological interest can be 
probed, in principle, through the detection of the redshifted, diffuse neutrino flux produced by 
all past supernovae — the so-called supernova relic neutrino (SRN) flux. In this work, we solve 
the SRN kinetic equations in the general case of three-generation flavor transitions followed by 
invisible (nonradiative) two-body decays. We then use the general solution to calculate observable 
SRN spectra in some representative decay scenarios. It is shown that, in the presence of decay, the 
SRN event rate can basically span the whole range below the current experimental upper bound — a 
range accessible to future experimental projects. Radiative SRN decays are also briefly discussed. 

PACS numbers: 14.60.Pq, 13.35.Hb, 97.60.Bw, 14.80.Mz 
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I. INTRODUCTION 



In the last few years, compelling evidence has emerged in favor of neutrino masses and 
mixings (see, e.g., [l| for recent reviews). Within a three-generation framework, the known 
flavor eigenstates of neutrinos and antineutrinos z/q, [a = e, /i, r) must then be considered 
as linear superpositions of mass eigenstates Ui {i = 1, 2, 3), 

i 

where the mixing matrix Uai (assumed here to be real for simplicity) can be expressed 
in terms of three mixing angles (^12,^13,^23) in the standard parametrization The z/j 
squared mass spectrum can be cast in the form 

^'+^'±^"^'1 ' (2) 

where and Am^ govern the two independent u oscillation frequencies, while fixes 

the absolute u mass scale. Current neutrino oscillation data (see, e.g., 0, Q) provide the 
best-fit values: 

W ^ 7.2 X 10"^ eV^ , (3) 

I ± Am^l ^ 2.0 X 10~^ eV^ , (4) 

sin^ ^12 ~ 0.29 , (5) 

sin^ ^23 ^ 0.50 , (6) 

and the upper bound (at 3a) 

sin^ ^13 < 0.067 . (7) 

The sign of Am^ is currently unknown, while is bounded from above by laboratory and 
astrophysical constraints (V < few x 10""^ eV 0). 

In general, massive neutrinos can not only mix, but also decay (see, e.g., the reviews in 
0, M, Q)- Neutrino decay has been invoked in the past, e.g., to solve the solar neutrino 
problem or the atmospheric neutrino anomaly (see [8| and references therein). These and 
other neutrino decay solutions have not been experimentally validated so far, implying that 
the decay lifetimes (r^ in the rest frame) must be sufficiently long to leave unperturbed 
the current neutrino phenomenology. Roughly speaking, this can be achieved by imposing 
TiE/niiL ^ 1 in the laboratory frame, where L and E are the neutrino pathlength and 
energy, respectively. Such arguments, properly refined and applied to solar neutrinos [L/ E 
O(IO^) km/MeV], lead to a safe lower bound 

— > 5 X 10"^ s/eV . (8) 
rrii 



^ This bound can be improved, for v ^ v decays, by one to two orders of magnitude (depending on specific 
scenarios) through the nonobservation of solar v^, v^. transitions in the KamLAND experiment, as 
recently reported in Ref. jlo| . 
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Similarly, from the supernova (SN) 1987A neutrino events [L/E ~ 0(10^^) km/Me V] one 
might naively guess a much stronger lower bound, 



>O(10^) s/eV, 



(9) 



which, however, does not hold for relatively large values of the mixing angle O12 as those 
implied by current data in Eq. (0) (see also the comment at the end of Sec. IV C). On the 
other hand, for neutrino decay effects to be observable in our universe, it must be roughly 
TiE/rrii ^ I/Hq, where Hq is the Hubble constant; by setting Hq = 70 kms~^Mpc~^ and 
taking E ~ 0(10) MeV (i.e., in the energy range probed by supernova neutrinos), a rough 
upper bound for the "neutrino decay observability" is obtained. 



The comparison of Eqs. (jHl) and ()10|) implies that many decades in Ti/rrii are still open to 
experimental and theoretical investigations. 

The exploration of relatively long neutrino lifetimes, and of the associated decay effects, 
can be effectively pursued through astrophysical and cosmological observations. For in- 
stance, evidence for neutrino decay might be signaled by an anomalous flavor composition 
of high-energy astrophysical neutrinos 125. Neutrino decay might also profoundly alter the 
flux of supernova relic neutrinos (SRN) and possibly push it close to the current experimen- 
tal SRN upper bounds [J^ J^], as recently pointed out in Ref. Given the promising 
prospects for detecting SRN signals in the Super-Kamiokande (SK) water Cherenkov detec- 
tor (if doped with Gd [3|) or in planned larger detectors such as the Underground Nucleon 
Decay and Neutrino Observatory (UNO) jl7| or Hyper-Kamiokande we think that the 
interesting idea of probing neutrino decay through SRN data y-5j appears as an opportunity 
which, although challenging, deserves further studies. 

On the basis of the previous motivations, in this work we aim at providing a general 
framework, as well as specific examples, for the calculation of nonradiative two-body decays 



(and for analogous antineutrino decays), where X is a (pseudo) scalar particle, e.g., a Majoron 
[7|. After considering the "standard" case of 3z/ flavor transitions without decay in Sec. II, 
we discuss in Sec. Ill the more general case of 3z/ transitions followed by decays, and give 
the explicit solution of the neutrino kinetic equations for generic decay parameters. Specific 
numerical examples (inspired by neutrino-Majoron decay models) are given in Sec. IV, in 
order to show representative SRN event rates and energy spectra in the presence of decay. 
Conclusions and prospects for future work are given in Sec. V. 

A final remark is in order. In this paper, we focus on invisible two-body decays [Eq. (Illj)] 
since they can be tested, in principle, by future SRN data. Invisible decays into three 
neutrinos are not considered here, also because the typical final-state neutrino energy would 
generally be too low to be probed by SRN observations. Radiative two-body neutrino 
decays, instead, do not suffer of such "low-energy" v detection problem, but are however 
severely constrained by the current astrophysical 7-ray phenomenology. A brief discussion 
of the photon background from hypothetical SRN radiative decays is given in Appendix A. 
Appendix B finally deals with some subdominant 3z/ oscillation effects in the atmospheric 
neutrino background, which are often overlooked in SRN searches. 



^ < 0(10^1) s/eV . 



(10) 



of SRN, 





3 



II. THREE-NEUTRINO FLAVOR TRANSITIONS WITHOUT DECAY 



In this Section we consider the standard case of flavor transitions among three stable 
neutrinos. Although the no-decay case for SRN has already been treated in previous works 
(see, e.g., [l^ 2^ 21, 2^, 13 recent contributions), we prefer to discuss its calculations in 
some detail, both to make the paper self-consistent, and to better appreciate the differences 
with the case of unstable neutrinos. In this Section we also define the notation, as well as 
our default choices and approximations for several inputs (from cosmology and astrophysics, 
neutrino oscillation data, and supernova simulations) needed to compute the SRN signal. 
The same inputs will be used in the presence of z/ decay (Sec. IV). 

We do not address here the problem of the uncertainties affecting these — and other equally 
admissible — input choices, being mainly interested to highlight the additional effects of neu- 
trino decay (as compared with the standard case). Careful considerations about central 
values and uncertainties of the various inputs will be necessary, however, when supernova 
relic neutrinos will be eventually detected by experiments and used to constrain models, 
with or without neutrino decay (see JisI] for a recent approach using simulated data). 

Finally, we remind that three-flavor oscillations affect not only the signal in SRN searches, 
but also the background. We discuss in Appendix B some subtle Su effects which are often 
neglected in calculating the atmospheric neutrino background rate. 



A. Input from cosmology and astrophysics 

The local (redshift z = 0) number density of mass eigenstates Ui per unit of energy, 
coming from all past core-collapse supernovae, is given by (see, e.g., ^19.]): 

n,XE)= dzH-\z)Rs^{z)Y,XE{l + z)) , (12) 

JO 

where R^^{z) is the supernova formation rate (per unit of time and of comoving volume 
at redshift z), Yy.{q) is the average yield^ of Pi for a typical supernova [per unit of initial, 
unredshifted energyg = E{1 + z)], and H{z) is the Hubble constant at redshift z which, in 
standard notation j^, reads 

H{z) = Ho[{l + zf{l + Qmz) - nAz{2 + z)]K (13) 

Since supernova relic neutrinos are ultrarelativistic (f ~ c), the above number density per 
unit of energy nj,^{E) can be identified, in natural units (c = 1), with the local relic z/j flux 
per unit of time, area, and eriergy. 

We assume, following jl^ Uj, that the function _Rsn is given by 



i?SN(^) =0.0122 i?sF(^)/MQ , (14) 



^ The yield of a species represents the integral of the v luminosity over the emission time, and is equal 
to the total u number times the normalized v energy spectrum. Since SRN signals are intrinsically time- 
averaged, we use V yields (rather than luminosities) throughout this work. Expressions for the 
functions are given below. 
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where the star formation rate Rsf is parametrized as 2^ 

RsA^) = 0.3 h,, . . Meyr-^Mpc-^ (15) 

45 + exp(3.8z) 

for z < 5 (and -Rsf = otherwi se Il9l| 'l. with /igs = Hq/65 kms^^ Mpc~^. (Other choices are 
also possible for Rsf] see, e.g., |22j.) 

The above expression for Rsf actually holds for an Einstein-de Sitter cosmology 
{flM, ^a) = (1, 0). For a different cosmology one has to apply a correction factor [l^ 

p R [(1+^)^(1 + ^^M^)-^^A^(2 + ^)]^ 

KsF -KsF -; -Ti , Uoj 

(1 + 2)2 

which cancels the dependence of riy- upon (^Mj^a) in Eq- f|T^- This cancellation does not 
apply in the presence of neutrino decay (see Sec. Ill B). When needed, we shall then fix 
VLm = 0.3, = 0.7, and h^^ = 1.077 (i.e., Hn = 70 kms~^ Mpc~^), consistently with recent 
determinations of cosmological parameters j25|. 



B. Input from neutrino oscillation phenomenology 

The smallness of the 5m^/Am^ ratio [see Eqs. © and (jl))] and of sin^ 6*13 [see Eq. (jTj)] 
lead to the approximate decoupling of the supernova Su dynamics into "high" {H) and 
"low" (L) 2z/ subsystems, governed by {kH,sm'^ 613) and (fc^, sin^ 6^12), respectively, where 
ku = -^^m? /2E and ki = 6m'^/2E are the two independent neutrino wavenumbers (see, 
e.g., 01 and references therein). Matter effects jl^ are governed by the potential V = 
±\^Gf Ne{x), where Ne{x) is the electron density at radius x (roughly falling as x~^ outside 
the supernova neutrinosphere, up to shock- wave effects jl^])' P^^^ (minus) sign refers 

to neutrinos (antineutrinos). The four cases corresponding to sign(Am^) = ±1 (normal 
or inverted hierarchy) and sign(V^) = ±1 (neutrinos or antineutrinos) lead, in general, to 
different physics. 

In supernova neutrino oscillations, it is customary to average out unobservable interfer- 
ence phases from the beginning, and to work directly in terms of level crossing probabilities 
Pij 2^1 . At the exit from the supernova, the yields of neutrino mass eigenstates Y^. are then 



given by 

yu. = T.p^^\Ua,\'yu. (17) 

(and similarly for antineutrinos), where ?7™- and Y^^ are, respectively, the mixing matrix 
elements in matter and the neutrino flavor yields at the neutrinosphere. 

The usual assumption Y^^ = Y^^ = Yp^ = Yp^ = Y^^ implies that the inner matrix product 
in Eq. (fT7|) depends only on Y^^, Yy^ for neutrinos (Fj,^, Y^^ for antineutrinos), and on the 
squared mixing matrix elements \U^\^. Such elements are given by 

i\UTi?. \U2?) = (cos^ On cos' OT^, cos' 6^, sm' ^^2, sin^ ^S) , (18) 



where 



cos2^r3 - ^""^^13 V/kH ^ -sign(V)sign(Am2) , (19) 



cos 2^13 - V/kHy + sin' 2613 

nnm COS 26*12 - V^/A;l . ,„„x 

cos26'"2 ^ , ~ -sign(y) , (20) 

(cos 2^12 - V/kiy + sin' 26^2 
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TABLE I: Relevant mixing matrix elements and supernova neutrino yields in the four possible 3i> scenarios. 
The k ^ I level crossing in the H system (if any) is reported in the 5th column. Sec the text for details. 



Scenario 


1 rrm|2 
\'~^el \ 


1 rrm|2 
\'-^e2\ 




k ^ 


/ 


Yield of i-th v state 














Y — Y 


v, normal hierarchy 








1 


3 ^ 


2 


Y,, = {1-Ph)Y,^+PhY,^ 

= PhY,^ + (1 - Ph)y,^ 














Y- — Y- 


normal hierarchy 


1 








— 




Yp2 = Y,y^ 

Ypa = Y,y^ 














Y — Y 


v, inverted hierarchy 





1 









Y — Y 














Y,, = (1 - Ph)Y,^ + PhY,^ 


u, inverted hierarchy 








1 


3 ^ 


1 


Y- — Y 

J^U2 — ^Vx 

= PhY,^ + (1 - Ph)y,^ 



and we have used the fact that \V/kH,L\ ^ 1 at the neutrinosphere. 

A further simphfication of Eq. (jl7p comes from the adiabaticity of the 1^2 level crossing 
in the L subsystem [P^ ~ for the mass-mixing values in Eqs. Q and (0)]- The matrix Pjj 
has then (at most) one nontrivial 2x2 submatrix with indices (fc, I) and off-diagonal entry 
Pff, where Ph is the k ^ I crossing probability in the H subsystem (equal for z/ and u (26l|). 
The final results for the yield of the i-th mass eigenstate at the exit from the supernova are 
collected in Table I. 

In general, the crossing probability Ph depends on the u potential profile V{x) (see, e.g., 
[2^). This dependence vanishes in the limiting cases of adiabatic transitions {Ph — 0) and 
strongly nonadiabatic transitions {Ph — 1), which correspond roughly to sin^ ^13 ^ 10~^ 
and sin^ 6*13 -C 10~^, respectively (see, e.g., |3fl|]). In our numerical examples, we shall 
consider for simplicity only the two limiting cases, ^ 

Pff = or 1 . (21) 

For our purposes, we can also neglect the secondary corrections due to Earth matter 
crossing [26|, |30| for arrival neutrino directions below the horizon. In this approximation, 
the number density of SRN in the flavor basis is simply given by 

n,^{E) = Y.\U^,fn,XE) , (22) 

i 

and similarly for antineutrinos. Finally, in the phenomenologically interesting case of super- 
nova relic z/e one gets, up to negligible terms of 0(sin^ 613), 

np^ {E) ~ cos^ n,, {E) + sin^ 6^2 n^2 (E) , (23) 

where sin^6'i2 is taken from Eq. (0). 



^ In general, Ph is a function of energy [20|, and possibly of time j2^. However, energy- and time-dependent 
effects on SRN signals are relatively small (as compared with the decay effects that we shall discuss) and 
are not considered in this work. 
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Supernova v yields (not redshifted) 
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FIG. 1: Supernova neutrino yields of Ve, Vg, and (curves labelled by Y^^, Yp^, Y^^, respectively), as 
defined through Eqs. (|^ ~ l(?7|) . 



C. Input from supernova simulations 

Concerning the three relevant flavor yields at the neutrinosphere = i^e,i^e,i^x), we 

assume for simplicity equipartition of the total binding energy Eb (taken equal to 3 x 10^^ 
erg) among all flavors, 

Y,{E) = -^^{E;{E),,a) , (24) 

where y:>{E) is the normalized neutrino spectrum (/ dE(p = 1), {E)^, is the average u energy, 
and a is a spectral parameter. A useful spectral parametrization is given in Ref. (3ll |: 



(a 



+ / E 



T{a + 1) \{E), 



~{a+l)E/(E)^ 



(25) 



where F is the Euler gamma function. From the admissible parameter ranges in Ref. j31|, 
we pick up the following values: 



((E),,,(E)p,,(E) 



,J = (12,15,18) MeV 



(26) 

a = 3 . (27) 

Figure 1 shows the functions Y^^{E)^ Yp^{E), and Y^^{E), which follow from the previous 
assumptions. We emphasize that such assumptions are not meant to provide "reference" or 
"best" neutrino spectra, but just a reasonable input for our numerical computations and for 
the comparison of results with and without decay. 
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FIG. 2: Supernova relic Ve spectrum (middle panel), and associated positron spectrum from Pe+P n+e~^ 
(right panel), in the presence of flavor transitions without decay. The spectra depend on the mass 
hierarchy (left panel), which can be either normal (NH) or inverted (IH). In the latter case, there is a 
further dependence on the crossing probability Ph- The representative cases Ph = ("large 6*13") and 
Ph = 1 ("small 6'i3") are shown. 



D. Results for supernova relic Uf. fluxes and positron spectra (no decay) 

Supernova relic z/g's can produce observable signals (positrons) through the reaction: 

Ue+P^n + e^ . (28) 

For our purposes, the cross section for this process can be approximated at zeroth order 
[3^ . with Ee+ = Ey — 1.293 MeV, where Ef,+ is the total (true) positron energy. We will 
generically refer to "water targets" but not to specific detectors; therefore, we will not 
include experiment-dependent details such as the difference between true and measured 
positron energy (i.e., the detector resolution function), or the detector efficiency function. 

Figure 2 shows our SRN fiux calculations for no decay, in both cases of normal hierarchy 
(NH) and inverted hierarchy (IH), as reported in the left panel. The middle panel refers 
to the absolute supernova relic Ve fiuxes [Eq. (j23p ] in units of i^e/MeV/cm^/s, as a function 
of neutrino energy. The right panel refers to the absolute positron event rate per MeV 
and per kton-year (kTy) in water, as a function of the true positron energy. In the case of 
normal hierarchy (solid curves), the Ve and positron spectra do not depend on the crossing 
probability Ph (s ee Table I). This case is indistinguishable from the case of inverted hierarchy 
with Ph = 1 j20|. In both such cases (NH with any Ph, or IH with Ph = 1) it is 

n,XE)= / dzH-\z)Rs^{z)[cos''e^2Y,XE{l + z))+siYi^e^2Y,SE{l + z))] . (29) 
Jo 

In the case of inverted hierarchy with Ph = (dashed curves in Fig. 2), the supernova relic 
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Ve spectrum depends only on Yy^^ 

n,XE)= dzH-\z)Rs^iz)Y,^{E{l + z)) , (30) 

^0 

and is thus peaked at somewhat higher energy [see Eq. and Fig. 1]. 

The results in Fig. 2 are in agreement with previous absolute estimates of the neutrino 
and positron spectra in water-Cherenkov detectors (see, e.g., [l^), modulo different choices 
for (some of) the inputs of SRN calculations. As we shall see, such results can be profoundly 
modified by neutrino decay. 



III. THREE-NEUTRINO FLAVOR TRANSITIONS AND DECAYS 

In this Section we discuss and solve the neutrino kinetic equations in the general case of 
Su flavor transitions plus decay. We consider generic rest-frame (anti)neutrino lifetimes Xj 
and associated decay widths,*^ 

r-' = r. = E r(z/. - y,) + T{u, z/,) . (31) 

Also, in this Section we do not set restrictions on branching ratios^ 

Bii^i Uj) = T{ui ^ Uj)/Ti , (32) 

and on neutrino decay energy spectra 

^,^^,^{E,, Ej) = PToh[ui{Ei) ^ Uj{Ej)] (33) 

(normalized to unity, / dEj ip{Ei, Ej) = 1). 

Notice that, for Ti/rrii values above the bound in Eq. (jSJ, supernova relic u flavor tran- 
sitions occur in matter (and become incoherent) well before neutrino decay losses become 
significant, so that hypothetical interference effects between the two phenomena |33] can be 
neglected. For our purposes, flavor transitions inside the supernova can thus be taken as 
decoupled from the subsequent (incoherent) propagation and decay of mass eigenstates in 
vacuum. 

An important remark is in order. In this work we consider only vacuum neutrino decays. 
It is possible, however, to construct models where fast invisible decays can be triggered 
by matter effects [3^ |3^ at the very high densities characterizing the supernova neutri- 
nosphere, even in the absence of vacuum decays (e.g, for diagonal neutrino-Majoron cou- 
plings; see, e.g. jl^). In such scenarios, matter-induced decays might thus occur before 
flavor transitions in supernovae, leading to a phenomenology rather different from the one 
considered in this paper (and subject to model-dependent constraints from supernova ener- 
getics 3^ 37, 38i,^, JO]). We emphasize that the results discussed in the following sections 



are generally applicable to vacuum neutrino decays occurring after flavor transitions, while 
possible matter-induced fast decays are beyond the scope of this work. 



^ In the laboratory frame, the widths are multipHed by the standard relativistic factor rrii/E. 
^ In the foUowing, we shall often use the symbol 'V' loosely, to indicate both neutrinos and antineutrinos. 
A distinction between ly and D will be made when needed to avoid ambiguities. 
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A. Kinetic equations 



In the presence of v decay, the local (2; = 0) SRN number density per unit energy n^.{E) 
depends on the history of all past {z > 0) decays having the z/j as initial or final mass 
state. The functions n^.{E, z) (number of Vi per unit of comoving volume and of energy at 
redshift z) can be obtained through a direct integration of the neutrino kinetic equations, 
as described below. 

The generic form of the kinetic equations for the phase-space distribution function / is 

m=C[f], (34) 

where C is the Liouville operator and C is the collision operator, embedding creation and 
destruction terms for the particle described by / (see, e.g., jMEIl)- 

For ultrarelativistic rehc neutrinos z/j (i? ~ p = | p^|), described at time t by riuXE.t) = 
Anp^ f B? {t) / R^, the Liouville operator takes the form^ 



The collision operator of unstable neutrinos reads 



n,,(E,t) . (35) 



where 



7D ■ 

C[n,XE,t)] = Rs^{t)Y,XE)+ E qAE,t) ~ n^^E^t) , (36) 

rrij >mi 

qj.{E,t)= dE'n,^{E',t)riu, u,) ^^,^^,XE',E) . (37) 

The right-hand side (r.h.s.) of Eq. contains two source terms and one sink term. 
The first source term quantifies the standard (decay- independent) emission of i>i from core- 
collapse supernovae. The second source term quantifies the population increase of Ui due to 
decays from heavier states; this term is absent for the heaviest neutrino state(s). The last 
(sink) term on the right-hand side (r.h.s.) of Eq. (jHUj) represents the loss of due to decay 
to lighter states with total width Fj; this term is absent for the lightest, stable neutrino 
state(s).'' 



B. General solution 

In order to solve the set of kinetic equations ()34j) - ()37p . we first rewrite them in terms of 
the redshift variable z = z{t) and of the redshifted energy e = e{E, z), 




Ro/R{t) - 1 

E/il + z) 



(38) 



^ R{t) is the universe scale factor for a Friedmann- Robertson- Walker metric, with H{t) = R{t)/R{t) in 
standard notation With respect to Ref. 0|, we drop a common factor {A'iTER^{t) / R^)^^ on both 



sides of the kinetic equations. 
^ Additional terms in the collision operator, due to Pauli blocking and to inverse reactions X -\- v 
can be safely neglected at the very low number densities of SRN. 



Vi 
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with associated partial derivatives 



dt ' 



\ dE J 



/ dz_ de \ 
I dt dt ^ 



dz de 
dE dE 



■de 



V 



1 d 
1+2 de 



(39) 



where the relation i = —(1 + z)H{z) has been used. With this change of variables, the 
Liouville operator depends on z only, and the kinetic equations can be directly integrated. 
More precisely, by defining the auxiliary function 



az)= r dz' H-\z'){l + z')-' 
Jo 



and the global source term 



Si{6{l + z),z) = Rst^{z)Y,Xe{l + z)) + ^ qj,{e{l + z), z) 

rrij >mi 

the neutrino kinetic equations ()34j) - (jnBjl can be cast in a compact form 



(40) 



(41) 



d_ r 

dz 



[l + z)n^^ e 



-miri^{z)/e 



SJe(l + z),z) 



which is easily integrated, giving 



1 



l + z 



dz' 



Hiz- 



-Si{e{l + z'),z')e-^^^^^^^''^-^^-'^^/' . 



(42) 



(43) 



By replacing back the variable E = e{l + z), one obtains the general solution of the neutrino 
kinetic equations. 



l + z 



dz' 



Rsn(z')Y,,^ IE 



1 



l + z' 



+ z 



l + z J 

-m,r,[i{z')^i{z)]{l+z)/E 



(44) 



which holds for generic neutrino mass spectra rrii, decay widths Fj, and decay energy spectra 
ipiy.^^^{E, E'). The effect of flavor transitions is embedded in the yields Y^^ (see Table I), 
in the same way as for the no-decay case discussed in Sec. II. Notice that the dependence 
of 1/H{z') upon the cosmology cancels with the factor Rsn{z') (see Sec. II A) but not with 
the other factors in Eq. (j44j) . Therefore, in the presence of decay, the SRN density acquires 
a dependence on the cosmological parameters {Qm,^a)- 

The double integration (over energy and redshift) implied by Eqs. ()37p and ()44j) can be 
performed numerically by following the decay sequence, i.e., starting from the heaviest state 
{qji = 0) and ending at the lightest state (Fj = 0). The observable local supernova relic Ui 
density is flnally obtained by setting 2; = in ny-{E^z). We conclude by noting that the 
case of no decay [Eq. (fT^] is obtained from the general solution [Eq. (jUj) ] at 2; = in the 
limit Tj — >■ 00 (i.e., F, = = g^j), as expected. 
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TABLE II: Branching ratios and energy spectra in 2v Majoron decay scenarios Vh ^ + X, for the 
extreme cases of quasi-degenerate (QD) and strongly hierarchal (SH) masses of the "heavy" and "hght" 
active neutrinos VhA- Analogous expressions hold for Ph decay, with the replacements Vk ^ Vk [k ~ h, I). 



Case 


Mass relations 


B{uh ui) 


B{vh vi) 


'tpvh-^ui {Eh, El) ipuh^ui {Eh, El) 


QD 


TJih ^ mi ^ mil — fTil 


1 





S{Eh -El) - 


SH 


mil S> ~ 


1/2 


1/2 


2Ei 2 / Ei\ 
'El E~hV~E~h) 



IV. APPLICATIONS TO SCENARIOS INSPIRED BY MAJORON MODELS 

In this section we apply the general results of Sec. Ill B to some representative decay 
scenarios, inspired by Majoron models. After a brief overview of the 2z/ decay case, we 
examine and compare a few representative Su decay cases, which provide SRN yields higher, 
comparable, or lower than for no decay. Simplificative assumptions will be made in all cases, 
in order to reduce the parameter space, and to highlight the main effects of neutrino decay. 



A. Two- family decay 

Let us consider a doublet of "heavy" and "light" neutrino mass eigenstates h'h,i {^k > ^i)- 
We briefly recall that nonradiative (invisible) decays of the kind 

Uh^'ui + X, (45) 

may arise through the coupling of Uh^i to a very light scalar or pseudoscalar particle X 
(assumed to be effectively massless for our purposes). In particular, X can be the Goldstone 
boson (Majoron) in models with spontaneous violation of the B—L symmetry of the standard 
electroweak model |4^ . 

There is a vast literature on neutrino-Majoron decay models (see, e.g., 0| and ref- 
erences therein) and related phenomenological constraints on neutrino-Majoron couplings 
(see, e.g., [I,i3£,^,^, 43i5i4^])- Since the neutrino-Majoron coupling can also contribute to 
the neutrino mass, model-dependent relations may arise between neutrino mass and decay 
parameters (see, e.g., [s^ 13). The branching ratios and final-state spectra for the two 
channels Uh vi and Vh — > are also functions of mo del- dependent couplings. 

Here, however, we do not commit ourselves to any specific theoretical model, and assume 
that the lifetime-to-mass ratio r/mh is a free parameter [subject only to the safe constraint 
in Eq. (jH)), when needed]. We also focus on two phenomenologically interesting cases in 
which the branching ratios and decay spectra become model-independent, namely, the case 
of quasidegenerate (QD) neutrino masses {ruh ^ mi ^ mh — rrii) and of strongly hierarchical 
(SH) neutrino masses {rUh mi ~ 0). Table II displays the relevant characteristics of the 
QD and SH cases (for ultrarelativistic neutrinos), obtained as appropriate hmits of the 
general results worked out in Ref. jist ]. 

In the QD case (Table II), the Uh decays only into z/;, which carries the whole Uh energy. 
When the decay is complete (i.e., for lifetimes rE/mh <^ Hq^ in the laboratory frame), the 
initial Uh energy spectrum is then directly transferred to the ui energy spectrum. 
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In the SH case, the Vh decays into either vi or vi with the same probabihty, but with 
different energy distributions (Table II). When the decay is complete, the initial Vh energy 
spectrum is then transferred to the final vi and vi spectra through convolutions with the 
'j/'-f unctions in Table II. For later purposes, in the SH case we define two (normalization- 
preserving) convolution operators V and acting over the Vh yield functions Yy^{E) as: 



Notice that, for supernova neutrino yields parameterized through Eq. |8ll |. the specific 
choice a = 3 in Eq. ()27|) makes the integrals in Eqs. pUj) and ()47|) analytical and elementary.^ 
Qualitatively, in the SH case the action of the operators T) and T) is to increase the neutrino 
yields at low energy. 

The above notation, although introduced for the 2v decay case, will be frequently used 
in the following subsections. In fact, we shall consider 2>y scenarios whose 2v sub-decays can 
be treated within either the QD or the SH approximation. 

B. Three-family decays for normal hierarchy and quasidegenerate masses 

In this section we consider a representative decay scenario which provides SRN densities 
generally higher than for no decay. The scenario involves normal mass hierarchy (+Am^), 
with masses much larger than their splittings [mf 3> Am^, 5m^). The quasidegenerate 
(QD) approximation, which forbids decays of neutrinos into antineutrinos and vice versa, 
can thus be applied (see Table II). In the context of supernova relic t'e's, we shall consider 
only antineutrino decays [u^ and z/2 — ^ z^i). 

For the sake of simplicity, we assume that -8(2/3 ^2) = B{us vi) = 1/2, and also 
that Ts/ms = T2/m2 = r/m. By construction, the decay scenario considered in this section 
is thus governed by just one free parameter^ (r/m). Notice that, for r/m ~ 0(10^") s/eV, 
SRN decay effects are expected to occur on a truly cosmological scale [see Eq. (fTUj) ]. For 
much larger values of r/m, the no-decay case is recovered. For much smaller values of r/m, 
SRN decay is instead complete, all SRN being in the lightest mass eigenstate Vi at the time 
of detection. 

Figure 3 shows the supernova rehc z/g energy spectrum, and the associated (observable) 
positron spectrum, for the decay scenario described above (and graphically shown in the left 
panel). The energy spectra for complete decay (red solid curves) appear to be a factor of 
~ 2 higher — and also slightly harder — than for no decay (black solid curve). This difference 
is entirely due to the role of z/^, as explained below. 

For complete decay, the final state is populated by stable i^i's only, coming both from 
the initial z/i component (Fp^ = YpJ and from fully decayed ^'2,3's {Yp^ + Ypg = 2Yy^, see 



For such choice, operators of the kind VD + VD and T)^ +1?^, which will appear in ?>v decay cases, also 
lead to analytical integrals (in terms of exponential integral functions). For generic, noninteger values of 
a, the analytical expressions become rather involved, and numerical evaluations are preferable. 
Oscillation parameters have been fixed previously, and the only relevant unknown (^13 or, equivalently, 
Ph) does not affect antineutrinos in normal hierarchy (see Table I). 




(47) 



(46) 
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V decay scenario energy spectrum e"^ energy spectrum 
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FIG. 3: (Color online.) Supernova relic v^. spectrum (middle panel), and associated positron spectrum 
from Ve ^ p ^ n + (right panel), for a decay scenario with normal hierarchy and quasidegenerate masses 
(left panel, with r/m and branching ratios assignments). The red and black solid curves correspond to the 
limiting cases of complete decay and no decay, respectively. The red dotted curves correspond to incomplete 
decay with r/m = 7 x 10^° s/eV. See the text for details. 



Table I), with unaltered neutrino energies (QD case in Table \\)}^ The final relic z/g density 
(given by cos^ ^12 times the final density of vi) is thus obtained by redshifting an initial 
yield equal to cos^ 6'i2lpg + 2 cos^ 6*12^^^. In the case of no decay for normal hierarchy, the 
initial yield is instead given by cos^ 612 Yp^ + sin^ OuYu:, [see Eq. fl^ ]. Therefore, while the 
Yp^ component is the same in the two cases, the weight of the Y^^ component for complete 
decay is 2 cos^ 9i2 — 1-42, much larger than for no decay, where the weight is sin^ 612 — 0.29. 
The stronger weight of Y^^ for complete decay leads to the increase in normalization, peak 
energy, and width, which is visible in Fig. 3.^^ 

For incomplete neutrino decay (i.e., for r/m ~ O(10^°) s/eV), one expects an intermediate 
situation leading to a SRN fiux moderately higher than for no decay. Figure 3 displays the 
results for a representative case (r/m = 7 x 10^° s/eV, red dotted curves), as obtained 
through the general solution of the kinetic equations worked out in Sec. III. Summarizing, 
the decay scenario examined in this section can lead to an increase of the SRN rate, as 
compared with the case of no decay. The enhancement can be as large as a factor ~ 2, the 
larger the more complete is the decay. 



The case of complete decay could also be obtained from the general solution in the limit r/m 
(derivation omitted) . 

The enhancement of the SRN yield for the decay scenario with normal hierarchy and quasidegenerate 
masses has also been discussed in Ref. llSl. 
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FIG. 4: (Color online.) Supernova relic Ve spectrum (middle panel), and associated positron spectrum 
from i>e + p — > n + e+ (right panel) , for a decay scenario with normal hierarchy and mi ~ (left panel, with 
r/m and branching ratios assignments). The green and black solid curves correspond to the limiting cases 
of complete decay and no decay, respectively. The green dotted curves correspond to incomplete decay with 
r/m = 7 X 10^" s/eV. See the text for details. 



C. Three-family decays for normal hierarchy and mi ~ 



In this section we consider a representative decay scenario which provides observable 
SRN densities generally comparable to the no-decay case. The scenario assumes that the 
mass hierarchy is normal (+Am^) and that the lightest state is basically massless (mi ~ 0), 
so that m2,3 ^ mi and the approximation of strong hierarchy (SH) can be applied to the 
decays of 1^2,3 (and of 1^2,3)- 

According to Table II (SH case), we can take B{v2 vi) = B(i'2 ^1) = 1/2- For the 
sake of simplicity, we extend such "branching ratio democracy" to all the z/3 decay channels, 
namely, we take B{v^ vi) = 5(^3 — > z/i) = -8(1^3 — >■ z/2) = B{i>3 — >■ 1^2) = 1/4. We also 
assume r^/m^ = T2/m2 = r/m, so that, as in the previous section, there is only one free 
parameter {r/m). The cases of no decay, incomplete decay, and complete decay, are then 
obtained for r/m much larger, comparable, or much smaller than O(10^°) s/eV, respectively. 

Figure 4 shows the supernova relic Ve and positron spectra for this scenario where, as 
depicted in the left panel, all decay channels are open. The results for the no-decay limit 
(solid curves) are identical to those in Fig. 3 and are not discussed again. The middle 
panel in Fig. 4 shows that the neutrino spectrum for complete decay (green solid curve) is 
significantly enhanced at low energy, as compared with the one for no decay. The positron 
spectrum, however, is very similar in the two cases (right panel). These features can be 
understood as follows. 

For complete decay, the final yield of stable z^i's comes both from the initial z/i's and from 
a complex decay chain, through the action of the operators T) and "D [Eqs. (I46p and (j47j) ]: 

vi^vi ^ ^1 , (48) 
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V2^yi^ \VY,, , (49) 

^ ^ Ivy,, , (50) 

1^3 ^ 1^1 ^ -VY,, , (51) 

p^^ui^ iw,3 , (52) 

V3^V2^i^i^ \VVY,^ , (53) 



\v 1 , , 

1^3 ^1^2 ^1^1 =^ ^VVY,, , (54) 

o 

U3^^2^iyi ^ l^VY,, , (55) 

o 

P3^i^2^i^i ^ -^VVY,, . (56) 

o 

This chain leads to a SRN density equation for complete decay formally similar to Eq. ()29p 
(no decay), but with the integrand in square brackets replaced by cos^ 6* times the sum of 
the terms in Eqs. (gH))-®: 

cos'ei2{Y^, + \vY^, + ... + \VVY,^ . (57) 



In the above expression for complete decay, terms which are linear and quadratic in the 
convolution operators produce a substantial enhancement of the SRN energy spectrum at 
low energy, visible as a "pile-up" of decayed neutrinos with degraded energy in the middle 
panel of Fig. 4.^^ At high energy, however, such terms happen to provide a contribution 
which is numerically close to the no-decay term sin^ 6'i2y^^ in Eq. ()29|1 . so that the high- 
energy tails of the SRN spectra for no decay and for complete decay (black and green solid 
curves in the middle panel of Fig. 4) are accidentally very similar to each other. Analogously, 
the case of incomplete decay (e.g., r/m = 7 x 10^° s/eV, green dotted curves), is appreciably 
different from the cases of no decay and of complete decay only at low energy. 

In this scenario, the interesting effects of decay are almost completely confined to low Ve 
energies, and are thus washed out in the observable spectrum, due to the cross section 
enhancement of high-energy features. In fact, the e"*" spectra for the three cases of complete, 
incomplete, and no decay, turn out to be very similar to each other (right panel of Fig. 4). 

In conclusion, from the results of this section we learn that there are neutrino decay 
scenarios which, despite a relatively complex structure (see Fig. 4), cannot be distinguished 
from the no decay case through future SRN observations, for any value of r/m above the 
safe bound in Eq. (jSJ. Similarly, such scenarios do not alter the "standard" (no-decay) SN 
1987A phenomenology, and thus provide a particularly clean example of how the naive limit 
in Eq. Q can actually be evaded. 



The terms linear and quadratic in T) and T) carry a mild dependence on the crossing probability Pa- 
through the neutrino yields YJ^^'s (see Table I for normal hierarchy). Figure 4 refers to the case Ph = 0; 
very similar results are obtained for Ph = 1 (not shown). 
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FIG. 5: (Color online.) Supernova relic z/g spectrum (middle panel), and associated positron spectrum 
from De + p n + (right panel), for a decay scenario with inverted hierarchy and generic ma (left panel, 
with t/tti and branching ratios assignments). The black and the blue curves correspond, respectively, to no 
decay and to incomplete decay with r/m = 7 x 10^" s/eV, for both Ph — 1 (solid) and Ph — (dashed). 
The horizontal (light blue) line at ~ corresponds to the case of complete decay. See the text for details. 

D. Three-family decay for inverted hierarchy 

We conclude our survey of Su decays by discussing a scenario where the SRN density 
is generally suppressed, as compared with the case of no decay. The scenario assumes an 
inverted mass hierarchy (— Am^), a fixed branching ratio B{u2 —>■ i^i) = 1/2 [= B{h'2 z/i)], 
and Ti/m = T2/m = r/m. Under these assumptions, only the decay U2 — > i^i (where the QD 
approximation is applicable) is relevant to SRN observations. In fact, decays to z/3 provide 
a negligible amount of UgS (oc sin^^ia). In particular, as graphically reported in the left 
panel of Fig. 5, the absolute value of makes no difference: It only opens (closes) the 
"dashed" decay channels for ~ (ms large), with no change on the relic flux for 
fixed B{h'2 i^i)- In the above scenario, the case of complete decay (r/m -C O(10^°) s/eV) 
is trivial: since the final state is populated only by z/3 (and 1/3), the relic density of is 
negligibly small (oc sin^^i3).^^ 

The nontrivial case of incomplete decay (r/m ~ O(10^°) s/eV) is then expected to 
lead to an intermediate suppression of the SRN density. Figure 5 shows the numerical 
results for the specific value r/m = 7 x 10^° s/eV, in both cases Ph = 1 (solid curves) 
and Ph = (dashed curves). The neutrino spectra for incomplete decay (blue curves) 
appear to be systematically lower than the corresponding no-decay spectra (black curves), 
although the difference is mitigated in the positron spectra (right panel of Fig. 5). A more 
substantial suppression of the positron spectrum can be obtained by lowering r/m (see next 
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This is true, in general, for complete decay in any IH scenario. 
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subsection). In conclusion, in the inverted hierarchy scenario of Fig. 5, the SRN signal is 
generally suppressed by neutrino decay, and can eventually disappear for complete decay. 



E. Overview and summary of 3u decay 

In the previous three subsections we have discussed, for three rather different scenarios, 
the cases of no decay [(r/m ^ O(10^°) s/eV], of complete decay [r/m <^ O(10^°) s/eV], 
and of incomplete decay for a specific value of r/m (= 7 x 10^° s/eV). We think it useful to 
show also the behavior of the SRN signal for continuous values of the free parameter r/m}^ 

Figure 6 shows the positron event rate integrated in the energy window Ee+ G [10, 20] 
MeV which, although lower than the current SK window might become accessible to 
future, low-background SRN searches For each scenario, the rate is normalized to 

the standard expectations for no decay and normal hierarchy (NH), and is plotted as a 
function of r/m. The black solid line represents the reference case (no decay in NH), which 
is indistinguishable from the case of no decay with inverted hierarchy (IH) and Ph = 1 (see 
Sec. II D). The black dashed line refers instead to the no-decay case with IH and Ph = 0. 
For no decay, variations in the hierarchy or in Ph appear to induce relatively small effect, 
which might be difficult to uncover once realistic experimental and theoretical uncertainties 
are considered in future SRN observations. Decay effects, however, enlarge dramatically 
the possible range of observable positron rates provided that r/m is in the cosmologically 
interesting range below 0(10^^) s/eV [See Eq. (fTUI) ]. For v decay in the quasidegenerate 
NH spectrum of Sec. IV B (red curve in Fig. 6), the positron rate rapidly increases for 
decreasing r/m, and reaches the "complete decay plateau" already at r/m ~ O(IO^) s/eV, 
with an asymptotic enhancement by a factor ~ 2.3. In the case of IH spectrum considered 
in Sec. IV D (solid and dashed blue curves in Fig. 6), conversely, the positron rate vanishes 
when approaching r/m ~ O(IO^) s/eV, for both cases Ph = 1 (solid) and Ph = (dashed). 
Finally, for the NH spectrum with mi ~ of Sec. IV C (green curve in Fig. 6), the positron 
rate appears to be almost the same as for no decay, at any value of r/m. 

Summarizing, neutrino decay can enlarge the reference no-decay predictions for observ- 
able positron rates by any factor / in the range ~ [0, 2.3].^^ 



Since the current experimental upper bound on the SRN flux from SK [13] is just a factor 
of ~ 2-3 above typical no-decay expectations [l^[l^ (including those in this work), future 
observations below such bound are likely to have an impact on neutrino decay models. If 
experimental and theoretical uncertainties can be kept smaller than a factor of two (a non- 
trivial task), one should eventually be able to rule out, at least, either the lowermost or the 
uppermost values in the range / G [0, 2.3], i.e., one of the extreme cases of "complete decay." 
Optimistically, one might then try to constrain specific decay models and lifetime-to-mass 
ratios through observations (a goal not reachable with current information ^3). Degeneracy 
between decay and no decay in specific models (as the one considered in Sec. IV C) will, 
however, set intrinsic limitations to SRN tests of neutrino lifetimes. 



This task implies extensive calculations through the general solution worked out in Sec. III. 
These results refer to a prospective positron energy window _Ee+ € [10, 20] MeV. For the current SK 
analysis window e [18,34] MeV, the range would be very similar, / e [0,2.7]. Variations of the 
reference inputs in Sec. II can also lead to wider ranges. 
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FIG. 6: (Color online.) Positron event rates in the energy range [10, 20] MeV for various decay scenarios, 
normalized to standard expectations for normal hierarchy and no decay. The r/m range in abscissa is well 
above the safe bound in Eq. 0. Notice how the expectations branch out (and then reach the complete 
decay limit) in the cosmologically relevant range r/m < 10^^ s/eV. See the text for details. 



V. CONCLUSIONS AND PROSPECTS 

Neutrino decays with cosmologically relevant neutrino lifetimes [rj/mj ^ 0(10^^) s/eV] 
can, in principle, be probed through observations of supernova relic z/g (SRN). We have 
shown how to incorporate the effects of both flavor transitions and decays in the calculation 
of the SRN density, by finding the general solution of the neutrino kinetic equations for 
generic two-body nonradiative decays. (Radiative decays are briefly commented upon in 
Appendix A.) We have then applied such solution to three representative decay scenarios 
which lead to an observable SRN density larger, comparable, or smaller than for no decay. 
In the presence of decay, the expected range of the SRN rate is significantly enlarged (from 
zero up to the current upper bound). Future SRN observations can thus be expected to 
constrain at least some extreme decay scenarios and, in general, to test the likelihood of 
specific decay models, as compared with the no-decay case. 
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In this work we have focused on theoretical SRN calculations, and have not attempted 
a phenomenological analysis of the sensitivity to neutrino decay through prospective SRN 
experimental data, as those which might be collected, e.g., in the proposed Gd-doped Super- 
Kamiokande detector in Hyper-Kamiokande or in the UNO project [17|. Such 
analysis would require a detailed characterization of the many uncertainties affecting both 
the background-subtracted SRN signal and the theoretical inputs, especially those related 
to supernova simulations and to the star formation rate. These uncertainties, although 
currently rather large, are likely to decrease in the future as more powerful supernova codes 
and new astrophysical observations will become available (including, hopefully, a galactic 
supernova explosion). When significant improvements will be made in this direction, a 
systematic approach to SRN predictions and to their uncertainties will certainly be required, 
in order to interpret future SRN data and to use them to constrain nonradiative neutrino 
lifetimes in a range [up to r/m ~ 0(10^^) s/eV] which is both largely unexplored and 
difficult to access by other means. 
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APPENDIX A: RADIATIVE DECAYS OF SUPERNOVA RELIC NEUTRINOS 



Radiative neutrino decays |46|, |47[ of the kind 



^i + l (Al) 

have been considered in a number of papers (see, e.g., 0, IEEE 3 references therein). 
Limits on the neutrino lifetime/mass ratio r/m for such decay modes can be derived from 
a variety of arguments 3|. E.g., one of the strongest bounds, 

r/m > 0(10'°) s/eV , (A2) 

is set by the nonobservation of infrared (IR) photons that would be emitted by Big Bang relic 
neutrinos [i^ and, in addition, by the nonobservation of scattering of high energy [O(TeV)] 
photons from distant sources on this hypothetical IR background [iP]. In the context of 
supernova neutrinos, the lack of excess 7 flux during the SN 1987A neutrino burst sets a 
limit r/m ^ 0(10^^) s/eV [44., 50]. In this Appendix, we estimate phenomenological limits 
on r/m from SRN radiative decays which, to our knowledge, have not been discussed so far 
in the literature. 

Hypothetical radiative decays of SRN produce a diffuse photon background, whose local 
{z = 0) density per unit of volume and energy is obtained by integrating over the (redshifted) 
source terms 



dz '"^^"^^ 



n,{E,) = l E qUE,{l + z),z) (A3) 
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where the gjj have an expression analogous to Eq. (|H7j). but with the appropriate final-state 
energy distribution for the photon in the laboratory frame, 



(for Ej > E^Am'jJm'j), where Am|j = m'^—mf, and the so-called asymmetry (or anisotropy) 
parameter a quantifies the amount of parity violation in the decay (3, EH {a E [— 1,+1] 
for Dirac neutrinos, while a = for Majorana neutrinos In calculating the photon 

density^^ from Eqs. ()A3|) and (|37p. one can assume at first order that the neutrino density 
at any z is basically "undecayed," namely 

„,(E,.).<(£,.) = ^f ^iSsKMn,(Bii|) . (A5) 

as derived from Eq. ()44|1 in the limit Tjmj/Ej 0. This assumption will be validated a 
posteriori below. 

In order to perform numerical calculations, the neutrino mass and decay parameters 
must be fixed. We assume a representative NH scenario with mi = 0, so that m|/Am^j ~ 1 
(strong hierarchy approximation) and the ip function in Eq. ()A4|) becomes the same for all 
decay channels. We also assume that T2/m2 = Ti/mi = r/m, so that the total photon fiux 
is 

rco rl^ rco ATT'' 

n,{E,) = {T/m)-'J^ -^J^^^^^iE',E,)^n'l^{E',z) , (A6) 

independently of the decay branching ratios (provided that they add up to unity, i.e., that 
only radiative decays occur). In the above equation, the inner sum extends to all the 
unstable states (both u and u). From Table I, it follows that the relevant neutrino yield to 
be integrated is 

+ ^3 + Y,, + = + , (A7) 

for any value of Ph- 

Figure 7 shows the results for the photon fiux n^{E^) in the above scenario, assuming 
r/m = 10^'^ s/eV and three representative values a = —1, 0, +1, corresponding to green, red, 
and blue curves, respectively. The curves cross at E^ ~ 1.21 MeV, where the a-dependent 
part of the integral in Eq. ()A6|) vanishes. Lower (higher) values of r/m would simply shift 
the curves upwards (downwards) along the y-axis. In the same figure, we superpose (as a 
dashed line) the best fit to the experimental 7 background fiux measured by the COMPTEL 
experiment for e [0.1, 30] MeV [sj, 

n'^^^'^iE.,) ~ 47r X (1.05 ± 0.2) x 10"' ( ) • ^^S) 

From Fig. 7 we derive that, in order to keep the 7 fiux from hypothetical SRN decay 
below the observed 7 background, it must be 

r/m > 0(10^") s/eV , (A9) 



For photons {v — c), the density (number of 7 per unit of volume and energy) also represents the flux 



(number of 7 per unit of area, time, and energy) in natural units. 
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7 flux from SRN radiative decay 




(MeV) 

FIG. 7: (Color online.) Energy spectrum of the photon flux coming from hypothetical radiative decays of 
SRN (solid curves), in a scenario with normal hierarchy, mi = 0, and T2/rn2 = T^/rris = 10^^ s/eV. The cal- 
culations are performed for three representative values of the decay asymmetry parameter a (= —1, 0, +1). 
The dashed line represents the power-law best fit to the 7 background measured by the COMPTEL exper- 
iment in the range e [0.1, 30] MeV 

up to an a-dependent factor of 0(1). Similar results would be obtained for inverted hier- 
archy and rris = (not shown). ^'^ Notice that this limit largely justifies a posteriori the 
"undecayed" approximation for n^, in Eq. flA5j) . 

We do not further elaborate the SRN constraint in Eq. ()A5|) . since it is not competitive 
with the one in Eq. ()A2|) . We note, however, that the supernova relic neutrino bound in 
Eq. ()A9|) is stronger than the one coming from the diffuse extragalactic stellar neutrino 
background [r/m > 0(10^^) s/eV] recently discussed in jHsj ]. 



Radiative decays of quasidegenerate neutrinos would instead produce 7's with typical energy much lower 
than in Fig. 7 for both NH and IH (not shown). We have checked that the corresponding bound on r/m 
would be degraded by roughly two orders of magnitude, as compared with Eq. ljA9|l . 



22 



APPENDIX B: REMARKS ON 2,v OSCILLATIONS OF THE ATMOSPHERIC 
NEUTRINO BACKGROUND 



In water-Cherenkov detectors, 3z/ flavor transitions affect not only the e"^ signal induced 
by supernova relic z/g, but also the background induced b y lo w-energy atmospheric neutri- 
nos. The irreducible components of the background [l3Lll4| are due to (1) interactions of 
atmospheric z/g and Vf.'-, and (2) decays of "invisible" /x^ (below the threshold for Cherenkov 
emission), induced by low-energy atmospheric and z/^. In Super-Kamiokande the 
background components (1) and (2) are characterized by parent neutrino energy spectra 
peaked at E ~ 100 and E ~ 150 MeV, respectively Q. 

In , the effects of neutrinos oscillations on the background are calculated by assum- 
ing pure 2u oscillations in the — >■ channel with maximal mixing (i.e., 5m? ~ 0, ~ 0, 
and sin^^23 — 1/2). Within the 2i/ approximation, the background component (1) is not 
affected by oscillations, while the component (2) is suppressed by the average effect of Am^- 
driven oscillations, through a factor = sin^ 2^23 sin^(l-27Am^L/£') ~ 0.5 sin^ 2^23 — 0.5, 
where [Am^] = eV^, [L] = m, and [E] = MeV. In terms of averaged flavor oscillation 
probabilities Pap, the 2i> approximation implies that {Pee, P^j,^, Pfj.e) — (1, 0.5, 0) [13^ . 

We remark, however, that the 2z/ approximation Sm^ ~ 0, usually valid for typical at- 
mospheric neutrino energies (^ 1 GeV), is not applicable in the SRN context. Indeed, for 
E rsj O(IOO) MeV and L ~ 0(i?®), the associated oscillation phase cannot be neglected 
[Sm^L/E ~ 0(1)]. Even for ^13 = 0, the low-energy atmospheric neutrino background is 
thus sensitive to the "solar" 5m^, rather than to the "atmospheric" Am^ (which can be 
effectively taken as 00). 

The possibility of probing (5m^-driven oscillations at very low atmospheric p energy 
(^ 1 GeV) would deserve, in itself, a separate investigation.^^ Here we only make some 
rough estimates of the relevant flavor oscillation probabilities Pa/B for 7^ 0. We assume 
isotropic atmospheric neutrino fluxes (produced at /i = 20 km from the ground level), equal 
components of u and u (with weights 3 and 1, respectively, to account for the different 
cross sections), and a u^/ue flavor ratio r = 2. Matter effects are estimated through a 
constant-density approximation for the Earth. The oscillations parameters are taken from 
Eqs. dni)-®, with either sin^^ig = or sin^ = 0.067 [see Eq. ((7j)]. We consider a repre- 
sentative neutrino energy E = 100 MeV, and average over all incoming u directions. 

Under these assumptions, and for sin^ 6^13 = 0, we find (Pee, P^i^i, P^le) — (0.77, 0.42, 0.11), 
which differ significantly from the values (1, 0.5, 0) obtained by setting 5m? = 0. The 
atmospheric electron events [background (1)] and invisible muon events [background (2)] are 
then modulated, respectively, by the factors fe — Pee + rP^e — 0.99 and ~ P^^ + Pee/f — 
0.48, which happen to differ only slightly^^ from the values (/e, f^) = (1, 0.5) considered 
in jl3l] by setting = 0. Analogously, for sin^ 6'i3 = 0.067 we find {Pee, P/i^i, P^ie) — 
(0.67, 0.17, 0.41) and {fe, fp) — (1.01, 0.50). Therefore, despite significant 5m^-induced 
effects on the oscillation probabilities, the irreducible atmospheric background rates do not 
appear to be significantly modified as compared with those obtained by (incorrectly) setting 
5m'^ = 0, independently of the value of sin^ ^13. 

In conclusion, the approximation 5m'^ = is not applicable, in principle, to the analysis 



For E > 0(1 ) GeV, subleading oscillations driven by Sm^ have instead been considered in several papers; 
see, e.g., and references therein. 

An accidental cancellation of effects is operative for the specific value r = 2. See, e.g., |5J|. 



23 



of the irreducible atmospheric neutrino background in SRN searches jlj]- However, in 
practice, this approximation does not appear to be harmful (according to our provisional 
estimates). More refined estimates of (5m^-driven oscillations might be of some interest in 
future high-statistics SRN searches 3) which are expected to find a signal above the 



irreducible atmospheric v background. These oscillation effects would instead loose interest 
if the "irreducible" background could be "reduced away" by tagging SRN events, as recently 
proposed in (l6l |. 
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